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We formulate the nonlinear isovector model in a curved background, and calculate the spherically 
- - - symmetric solutions for weak and strong coupling regimes. The usual belief that gravity does 

not have appreciable effects on the structure of solitons will be examined, in the framework of the 
, calculated solutions, by comparing the flat-space and curved-space solutions. It turns out that in the 

^ ■ strong coupling regime, gravity has essential effects on the solutions. Masses of the self-gravitating 

solitons are calculated numerically using the Tolman expression, and its behavior as a function of 
^ |. the coupling constant of the model is studied. 
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I. INTRODUCTION 



There has been considerable interest in the localized solutions of the Einstein's equations with nonlinear field 
sources in recentyears (0, 0, ||]). Gravitating non-abelian solitons and black holes with Yang-Mills fields is 

investigated in [|6|. Such problems were not investigated earlier in the history of GR, mainly because of two reasons: 
' 1. It was widely accepted that the gravitational effects are too weak to affect -in an essential way- the properties of 
, soliton solutions of nonlinear field theories. 2. The resulting equations are usually formidable such that the ordinary 
0^ ■ analytical approaches become idle. More recently, however, the availability of high speed computers and advanced 
numerical methods have changed the case, and extensive numerical attempts have been made in this direction (see 
• e.g. 387N Term Project in It has emerged from recent studies that the effects due to the inclusion of gravity 
■ are not always negligible. Consider, for example, the Einstein- Yang-Mills (EYM) system. It has been shown that the 
O EYM equations have both soliton and black hole solutions (|l|l, and ||]). This is in contrast to the fact that vacuum 
Einstein and pure Yang-Mills equations do not have by themselves soliton solutions. We can therefore conclude that 
$H I gravity may have dramatic effects on the existence or non-existence of soliton solutions of nonlinear field equations. 

Another interesting example is the discovery that black hole solutions may have Skyrmion hair . It was previously 
^ \ believed that stationary black holes can only have global charges given by surface integrals at spatial infinity (the 
so-called no- hair theorem). 

. _ In the ordinary 0(3) model, spherically symmetric solutions have an energy density which behave like l/r^ at large 
^ , distances (0). When formulated in a curved background, this model leads to a spacetime which is not asymptotically 
fiat, and the ADM mass is not well defined ((|]). A nonlinear 0(3) model (thereafter referred to as the isovector 
model) was introduced in ([|lO|), which possesses spherical, soliton-like solutions with a energy density at large 
distances. Such a model, is therefore expected to be well behaved in an asymptotically fiat background. In the present 
paper, we examine this model, and discuss its self-gravitating solutions. These new solutions are compared with those 
obtained previously in a flat spacetime. 

The present manuscript is organized in the following way. In section ||, we will review the isovector model of . 
In sectionpR flat-space solitons of the isovector model and their resemblence to charged particles are introduced. In 
section |IV|, the isovector model will be reformulated in a curved background. The resulting differential equations for 
a spherically symmetric ansatz will be introduced in this section, together with the necessary boundary conditions. 
These equations will be solved numerically, for several choices of the coupling constant. We will compare the self 
gravitating solutions with those obtained for a fiat spacetime. Soliton masses using the Tolman formalism will be 
discussed in section ^ together with the behavior as a function of the model parameter. Section VI will contain the 
summary and conclusion. 
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II. ISOVECTOR MODEL 

Consider an isovector field (pa (o — 1, 2, 3) with a 5^ vacuum at 



(1) 



Each component (j)a is a pseudo-scalar under spacetime transformations, and 0o is a constant. A topological current 
can be defined for such a field according to (pX[|) 

= ^e''''"''eabcd,^ada<l>bdp^c- (2) 

For the time being, spacetime is assumed to be the flat Minkowski spacetime and fj,, ... — 0, 1, 2, 3 with x° = t {c = 1 
is assumed throught this paper), e^'^"'^ and eabc are the totally anti-symmetric tensor densities in 4 and 3 dimensions, 
respectively. It can be easily shown that the current (^) is identically conserved {d^Jf^ = 0), and the total charge is 
quantized 

rd^x ^ ± / ^dS, = ne, (3) 
27r / dSx 

where e = 20^. In this equation, dSx and dScf, are area elements of 5^ surfaces in the x-space (as r — > oo) and (/)-space 
(as (?!) — > <^o), respectively. The current (^) can identically be written as the covariant divergence of an anti-symmetric, 
second-rank tensor 

d^F^"" = r, (4) 

where 

p^u ^ [e^^^^^dMpcl^c + dfiC^] , (5) 

ZTT 

in which Ba is an auxiliary vector field. The dual field * F with the tensorial components 

*F^'^ = ]-e^^''''PF^p = ^(2e„,,^, 9^-0,9-^, + d'^C^ - d^n, (6) 
Z 4:7r 

satisfies the equation 

*F>''' = 0, (7) 
provided that the vector field is a solution of the following wave equation 

□C^ - a''(a„C") = 2eabcdu{(i^ad''cj)bd^'cl^,). (8) 
It can be easily shown that the right hand side of this equation defines another conserved current 

= 2eabcdc.{cbad''<f>bd"cf>c), (9) 

with 

d^K^ = 0. (10) 
Using the language of differential forms, (^ can be written in the following form 

F = G + H, (11) 
where the components of the 2-forms G and H are given by 

Gfj.u = -^e"^ ^^eabc<Pada(l}bdi3(l}c, (12) 



and 



H^. = ^e^.^f'dpB^. (13) 
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We now have 

dF = 0, (14) 

and 

d*H = 0. (15) 

The 2-form F is therefore Hodge-decomposable, and cohomologous with G (i.e. they belong to the same cohomology 
class, since they differ only by an exact form). The resemblance of equations (|) and (^ to the Maxwell's equations 
and the capability of this model to provide non-singular solutions behaving like charged particles were discussed in 
pO|. In the next section, we will only outline the main results valid in a flat spacetime. 



III. FLAT SPACE SOLITONS 

The requirement of having non-singular, finite energy and stable solitons, severely restrict the possible choices of 
the lagrangian density of the isovector field. Let us follow ||l^, and adopt the following lagrangian density which 
satisfies the above requirements: 

C^Xid^^ad'^^af -bo{l-^f, (16) 

(Po 

with A < 0, and bo real constants. The potential V{(f>) = bo{l — <f>/(l)o)'^ satisfies the following conditions 



5^ 



^(^o) = 0, ( — )^„ = 0, and (Trr^)0„ > 0, (17) 



which leads to the spontaneous breaking of the (global) 50(3) symmetry of the system. The dynamical equation for 
the isovector field is easily obtained, using the variational principle 6 J Ld'^x = 0, which leads to 



d^dMd''cj)ad''cj)i + d^d^d^cj^bd'^cj^a + d,cl>bd''^bOcf,a = - ■ (18) 



Similar to the ansatz used in the Skyrme model ([Q)' start with the so-called hedgehog ansatz 

0a = 0(r)-, (19) 
r 

where x", a = 1, 2, 3 represent the Euclidean coordinates x, y, and z, respectively. This ansatz immediately leads to 

Ki" = Q, = Q, B = 0, (20) 

and 

where f is the unit vector in the radial direction, and e is the elementary topological charge defined in (^. Note that 
F°^ — Ei and F'^ = —CijkBk- The charge density corresponding to this ansatz is easily obtained to be 

JO 3 02 

27r dr 

which leads to 

' pAnr^dr = 20j^ = e, (23) 

showing that the ansatz bears unit topological charge. Note that in deriving this result, we have used the boundary 
conditions 

Mr = 0) = 0, (24) 
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and 

4> ^ (fio, as r — > oo. 
It can be shown that the foUowing asymptotic solutions are vahd: 



r 1 
ao 



+ 



close to the center of the soliton (r 0), where ai is a dimensionless constant, and 

(j){r) ~ (f>o 

and 



(25) 

(26) 

(27) 
(28) 



far from the soliton (i.e. r oo). 

It can be seen that the total soliton energy 



M = / T°d^x = 



Anr^dr = 5.210^ [bo{-4:X)^] 



1/4 



(29) 



Using the scale transformation r —f ar (while keeping (j) unchanged), it can be shown that M{a) has a minimum at 
a = 1, which is a signature of the stability of the soliton under radial perturbations. 



IV. SELF-GRAVITATING ISOVECTOR SOLITONS 



By self- gravitating isovector solitons, we mean static solutions of the coupled isovector-gravitational equations which 
are everywhere regular and represent localized lumps of energy. By numerically integrating the coupled nonlinear 
equations, we will show that such solutions do arise depending on the value of the model parameters. Based on or 
results, we will also criticize the widely expressed view that gravity has only a minute effect on the structure and 
properties of extended solitons. 

Let us start with the action 



A 



R 



IQttG 



C 



M 



-gd X, 



(30) 



in which G is the gravitational constant, R is the curvature scalar. Cm is the lagrangian density of the matter source, 
and g is the determinant of the metric tensor. As the source of gravity, we consider the isovector field ([l6|). By varying 
the action (pO|) with respect to and 0a, we obtain the corresponding field equations: 



Rfiiy — SttGA 



V{cl>) 



and 



{d^c^bd0c^b)o<t>a + g''''d,,^adAd^(t>bdp^b) 



1 <l>a dV{^) 

4X (j) d(j) ■ 

By contracting equation (|3ll) , we obtain the following equation for the curvature scalar: 

R ^ -327rG\/(0). 



(31) 



(32) 



(33) 



This equation is useful, since it expresses a simple relation between the curvature scalar and the self-interaction 
potential of the isovector field, and provides a means to check some of the calculations. We employ the coordinates 
= {t,r,9,(t)), and the general, spherically symmetric, static metric = diag(— B(r), A(r), r^, r^ sin6'). For the 
hedgehog ansatz (l9h, the independent field equations become 



B A \ 2A r 



(34) 
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B' A' r^.frc^" Acj,^ rAV{^\ 2A 2 

— = -647rGA — ^ H — H 

B A \ AA AX J 



(35) 

r r 



3(t)'^(j)" 2(f>^(j)" 3A'(j)'^ B'(j)'^ 2(/)'3 200'2 A'<j)'^(j)' B'cp^cj)' Acj)^ _ 1 dV 



and 



Note that equation (|36| ) reduces to the flat-space equation (|18D, by putting A — B — I. Equations ( ^41 ) and (|3 
however, become inconsistent for obvious reasons (the spacetime cannot be flat in the presence of matter sources) 
Let us introduce the following three length scales 



-4_A^ 
bo 



(37) 



1 



rgi = -r^===, (38) 



and 



rg2 = (39) 



which appear naturally in equations ( |34| ) to (|3q). It can be seen that Tq is proportional to the geometric mean of rgi 
and rg2 : 



(40) 

Emergence of new length scales is similar to what happens in non-abelian gauge fields (Q]), and leads to the 
appearance of new branches of spherical, static solitons. Having additional length scales, one expects departures from 
the Einstein-Maxwell solutions and the appearance of more subtle features. Using these two length scales, equations 
(^) to (|36|) can be made dimensionless and suitable for numerical integration; 



B' A' 1 /xu 



B A e \ 2A 



,'4 



(41) 



B' A' 1 /xu"^ Au^ \ 2A 2 

— -xAW{u)]+ , (42) 



B A e\ AA 



X X 



and 



Zu'^u" 2u^u" ZA'u'^ B'u'^ 2u'^ 2uu'^ A'u^u' B'u^u' _ dWju) 

A^ ^ 2A3 ^ 2A^B ^ ^ ^ x^A ^ x^A^ ^ x^AB du ' 

In these equations, u — 4'/4'o, x — r/vo, and t=a^ which is the dimensionless parameter of the system: 



(43) 



(44) 



To J 32TTG(j)lV^Xbo 

and 

W{u)^{l-uf. (45) 
The asymptotic behavior of the non-singular solutions can be found by using the following series ansatze: 

oo oo oo 

u{x)^^Cix\ A{x)^^aix\ B{x)^^bix\ (46) 

i=0 1=0 i=0 
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r ^0 


r — > OO 


ao = 1 

1 9uj*-4 
«2 24 f 
1 " 

«3 = To7 

1 12555t^l''-12240cj'^ -2880cj*e+2000iJ^ -936f 


«2 = -^ 

/ 1 

"4=4^^ 

„/ 1 5-16e^ 
«6 40 e3 
„/ 5-326^ 


"4 72000 /^i.:''! 


"8 80 


bo = 1 

bo - 1 ^'^''+4 
"2 24 e 

I, 13 w 

"3 - -307 

I 1 1035u^'' + 220u'' + 1440a;*'e+268e 
"4 8000 


U _ 68 
"10 - 45e 
U _ 96 
0l2 - 5572- 


Cl = CiJ 

C2 = -Tob^ 

1 45tjl"-35ij''+30a;*'e-9e 
•^3 — 750 ijE>£ 

_ 1 35505a;l"+6260u'' + 10920a;*'e-3456e 
''4 1350000 etj» 


4 = -2 
C8 = -28 

„' _ 28 
ClO — T 

c'i2 = -1648 



TABLE I: Leading, non- vanishing coefficients of the asymptotic solutions. 
O TTl I 




o 



1 2 3 4 5 6 

X 

FIG. 1: Variations of the isovector field amplitude for several values of the parameter e. Downward: flat space, e=0.2,0.1, and 
0.05. 

for X — > 0, and 

OC / OO f OO J/ 

-(-) = i+E^' ^(-) = i+E|' ^(-)-i + E^' (47) 

i—1 i—1 i—1 

for X 00. The unknown coefficients are calculated by inserting these ansatze into equations ( ^ ) to (^3|) and 
balancing the terms of the same order in x. Table | shows the leading coefficients for the above asymptotics. 

In order to solve the coupled nonlinear DEs numerically, we used the Gerald's shooting method, which is based on 
two guesses for the initial slope of the unknown functions. Using these initial guesses, the equations are integrated 
via the Runge-Kutta-Fehlberg method, reaching the end point of the independent variable. A better guess for the 
initial slope is then found by comparing the end point values with the boundary conditions, and interpolating for the 
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FIG. 2: Variations of the metric coefficient A for several values of the parameter e (=0.5,0.2,0.1, and 0.05 downward). 

initial slope ([p^). The boundary conditions for asymptotically flat spacetime read 

u^O as a; 0, (48) 

and 

u ^ 1, A ^ 1, and _B — > 1, as x ^ oo. (49) 

The procedure described above is iterated until the correct boundary values are reached with a reasonable accuracy. 

In order to test the method, flat space solutions were flrst computed and compared to the solutions obtained via 
energy minimization algorithms ( ||l0[| ). Figure ^ shows numerical variations of the u{x) function for several values of 
the parameter e. The corresponding results for the metric coefficients A{x) and B{x) are shown in Figures |^ and |^, 
respectively. 

It is seen that for e of the order of unity, the self-gravitating solutions differ only slightly from the flat-space solution. 
The difference vanishes completely as e — > oo. As e becomes much smaller than 1, significant differences with the 
fiat-space solution emerge. For example, the metric signature changes at some r for sufficiently small e. 



SOLITON MASSES 



Although there are still controversies about an exact definition for the total mass of a self-gravitating system ( [g5| , 
16[,|1^), we adopt the Tolman formalism for computing the total mass of the gravitating solitons; 



Mt = J IrV^d-'x, (50) 

where 

It - - - r| - T|, (51) 



Z5 




1 2 3 4 5 6 

X 

FIG. 3: Variations of the metric coefiicient B for several values of the parameter e (=0.2,0.1,0.5, and 0.05 downward). 



and 




-e — , , — () 



10 10 



10 10 10 

e 



10 10 



FIG. 4: Variations of the integral St as a function of the dimensionless parameter e. 
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(52) 



In most cases (including the model we are discussing here), the second term is absent since the matter lagrangian 
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does not contain derivatives of the metric tensor. For the isovector model, we have 



It = 4A 



2(/)2\ /2r0'0" rA'(j)- 



A 



A 



'^^^ ^ :,9V (rA! rB' ^ 
~2r(t>' — -[ — + — +6 



(53) 



Using the transformations ( [44| ) and ( |45| ) , the total mass can be written in terms of a dimensionless integral: 

MT = 7T^l[bo{-^Xf]'/^ST, (54) 

where 

/■°° 2^2 2^2 \ /2xu'w" xA'u'"^ u'^ Auu' 4^^ 
•^^--4/ + — + ______ 



A 



A 



X X- 



xA' 



xB' 



W{ti) 



u 
~A 



2u^ 



2xu' ' )VABx^dx. 
ou 



(55) 



and rj = — 4A. The dimensionless integral St were computed numerically for several values of the parameter e. The 
results are shown in Figure |[ 

It is seen that (1) St — > 1-66 as e ^ oo, leading to the asymptotic (flat space) mass, (2) There is a maximum mass 
around e ~ 1, and (3) The total mass decreases as e ^ 0, with Mt ~ at e ~ 0.055. It is also interesting to note 
that in the asymptotic series solutions summarized in Table |[ the coeflticients and 6q vanish, which imply vanishing 
total mass of the soliton as deduced from the asymptotic form of the metric [|l9| . 



VI. SUMMARY AND CONCLUSION 



We extended the isovector model to incorporate the effects of gravity. The resulting equations were integrated 
numerically for spherically symmetric ansatz, using the Gerald's shooting method. It was found that for large values 
of the dimensionless parameter of the system the effect of gravity is negligible. For small values of e, however, gravity 
has a considerable effect on the qualitative and quantitative behavior of the solutions. Such dramatic changes in the 
behavior of the spherical solutions in the presence of gravity were also reported in the framework of EYM equations 
(0, 1^, and [^). Gravitating solitons of the isovector model in an asymptotically flat background bear quantized 
topological charges, exactly similar to the flat-space solitons. The quantization is due to a tt^ homotopy between 
the boundary of the curved space [S'^ at r — > oo), and the vacuum S"^ of the isovector field [(j)a4'a — 0o)- This is 
in analogy with the quantization of the magnetic pole intensity in the t'Hooft Polyakov monopoles ([|l8j). Solutions 
presented in this paper (Figures ^ to |^) do not exhibit horizons. Using the well-known result from general relativity 
([0), T2/T1 = ^B{r2)/B{ri), where n and T2 are the emitted (at ri) and detected (at 7-2) periods of photons, 
event horizons correspond to B{ri) = 0, which are not fulfilled by the present solutions. However, the appearance of 
horizons and signature changes is not ruled out and will be addressed elsewhere. In particular, it should be interesting 
to know whether isovector black holes have hair. As pointed out in solitons can make bound states with ordinary 
black holes to form hairy black holes. In such a case, the total ADM mass of the hairy black hole is the sum of the 
mass of the bare black hole, the mass of the soliton, and the gravitational binding energy between the two (||l3||). 
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